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Abstract 

The Nonlinear Schrodinger (NLS) equation is used to model surface waves 
in wave tanks of hydrodynamic laboratories. Analysis of the linearized NLS 
equation shows that its harmonic solutions with a small amplitude modulation 
have a tendency to grow exponentially due to the so-called Bcnjamin-Fcir in- 
stability. To investigate this growth in detail, we relate the linearized solution 
of the NLS equation to a fully nonlinear, exact solution, called soliton on finite 
background. As a result, we find that in the range of instability the maximum 
amplitude increase is finite and can be at most three times the initial amplitude. 

Keywords : Nonlinear Schrodinger equation, Benjamin-Feir instability, 
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1 Introduction 

This is an initial work on 'Extreme Waves in Hydrodynamics Laboratory'. Extreme 
waves here refer to very high amplitude, steep, waves that can appear suddenly from 
a relatively calm sea. Although extreme waves are very rare and unpredictable, they 
are still very dangerous to ships in case they meet. 

We model the problem of extreme waves using dispersive wave modes. The spe- 
cific property of dispersion is that waves with different wave length propagate with 
different phase velocities. In the following we assume that the wave field has a fre- 
quency spectrum that is localized around one frequency. Then the envelope of the 
wave field is described by the NLS equation [3] . The NLS equation is an amplitude 
equation for describing the change of envelope of a wave group. This equation is very 
instrumental in understanding various nonlinear wave phenomena : it arises in stud- 
ies of unidirectional propagation of wave packets in a energy conserving dispersive 
medium at the lowest order of nonlinearity [3]. 
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In this paper we analyze the behavior of the wave group envelope using both 
linear and nonlinear theories. Linear theory predicts exponential growth of the 
amplitude when certain conditions are satisfied — the Benjamin-Feir instability [1]. 
However, when the amplitude becomes large nonlinear effects must be taken into 
account, that, as it turns out, will prevent further exponential growth. The aim of 
this paper is to find the maximum amplitude of waves when the amplitude growth is 
triggered by the Benjamin-Feir instability. Also, we investigate how this maximum 
amplitude depends on the growth rate parameter from the Benjamin-Feir instability. 

2 Modelling of Waves Envelope 
2.1 Linear Theory 

In linear theory of water waves, we can restrict the analysis of a surface elevation 
rj(x, t) to one-mode solution of the form rj(x, t) = a e 4 ( fca '~ wt ) + c.c. (complex conju- 
gate), where a is a constant amplitude, k is wavenumber, and uj is frequency. Then 
a general solution is a superposition of one- mode solutions. For surface waves on 
water of constant depth h, the parameters k and uj are related by the following linear 
dispersion relation [3]: 

uj 2 = g k tanh k h, (1) 
where g is gravitational acceleration. We also write the linear dispersion relation as 
uj = £l(k) = k W 1 9 tai ^ h ^ A, This dispersion relation can be derived from the linearized 
equations of the full set of equations for water waves. The phase velocity is defined 
as and the group velocity is defined as ^ = fi'(fc). 

Since the linear wave system has elementary solutions of the form g'f^-^W^ 
it is often convenient to write the general solution of an initial value problem as an 
integral of its Fourier components [3] : 

/oo 
a {k) e i{kx-mt) dk ^ ( 2 ) 
oo 

where a(k) is the Fourier transform of rj(x, 0). Writing dispersion relation Q(k) as a 
power series (Taylor expansion) about a fixed wavenumber ko and neglecting 0{k?) 
terms, we find 

(3 k 2 = n(k + as) - n(k ) - n'(k ) k, (3) 

where /3 = \ fl (ko). Let us define k = ko + k, t = t, and £ = x — £l'(ko) t. Then 
equation ^ can be written like 

/oo 
aiko + ^e^^e^^^dK. (4) 
oo 

Denoting the integral in Q with i/j(£,t), we find that 4 ) (£,i T ) satisfies 



Investigation of the maximum amplitude increase 



3 



This is the linear Schrddinger equation for narrow-banded spectra. Equation ([5]) 
is a partial differential equation that describes time evolution of the envelope of a 
linear wave packet [4]. Equation ([5]) has a monochromatic mode solution ip(£,i T ) = 



2.2 Nonlinear Theory 

Assuming narrow-banded spectra, we consider the wave elevation in the following 
form r}(x,t) = i/t(£, r) e <(*o*-"o*) + c .c, where r = t, £ = x - fi'(Jfeo)t, is 
a complex valued function (called the complex amplitude), k^ and wo are central 
wavenumber and frequency, respectively. The evolution of the wave elevation is a 
weakly nonlinear deformation of a nearly harmonic wave with the fixed wavenumber 
ko. If we substitute rj to the equations which describes the physical motion of the 
water waves (see below), then one finds that the complex amplitude ip(£, r) satisfies 
the nonlinear Schrddinger (NLS) equation. 

As example, the NLS equation can be derived from the modified KdV equation 
n t + iQ (—id x )n + | d x rf = [2]. With r = t and £ = x — O'(fco) t, the corresponding 
NLS equation reads 

+ M 2 V> = o, ft 7 effi, (6) 



where /3 and 7 depend only on k$ : 

P = \ri{**\ (7) 

= _9 / 1 fc 
7 4 VO'(*o) - ^'(O) 2fi(A;o) - fi(2fc ) 
where w = fi(fe) is the linear dispersion relation [6]. The coefficients j3 and 7 in this 
paper have opposite signs compared to the corresponding coefficients in [6]. This 
equation arose as a model for packets of waves on deep water. 
There are two types of NLS equations : 

• If /3 and 7 have the same sign, i.e. /3 7 > 0, then ^ is called the focusing NLS 
equation (an attractive nonlinearity, modulationally unstable [Benjamin-Feir 
instability]) [5]. 

• If /3 and 7 have different sign, i.e. ^7 < 0, then ^ is called the de-focusing 
NLS equation (a repulsive nonlinearity, stable solution) [5]. 

The wavenumber fc cr it, for which j3 7 = holds, is called the critical wavenumber 
or the Davey-Stewartson value. Using the physical quantity g = 9.8 m/s 2 and the 
water depth h = 5 m, then for wavenumbers ko > k CT i t = 0.23, the product ^7 > 0, 
and the NLS equation is of focusing type. The corresponding critical wavelength is 
27.41 m. 
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In the following, we consider only focusing NLS equations. The NLS equation 
^ has a plane-wave solution 



„2 . 



A(t) =r e^ r o T . (9) 
In physical variables, the surface wave elevation is given as rj(x,t) = 2rocos (kox — 

1-7 
ko 



uj^t + jr^t). Note that the corresponding phase velocity is fc J r ° . In the following 
section, we analyze the stability of this plane-wave solution. 



3 Benjamin— Feir Instability 

To investigate the instability of the NLS plane-wave solution, we perturb the £— 
independent function A{t) with a small perturbation of the form e(£, r) = A(r) B(£, r). 
We look for the cases where under a small perturbation the amplitude of the plane 
wave solution grows in time [3]: 

^,t) = A(t)[1 + B^,t)}. (10) 

Substituting (10) into ([6]), and ignoring nonlinear terms we obtain the linearized 
NLS equation : 

iB T + pB^ + 7 rl(B + B*) = 0. (11) 
We seek the solutions of ( |11[ ) in the form 

B(£, t) = B ie ^ T+iK ^ + B 2 e( a * T - iK t\ (12) 

where Bi,B 2 E C, ko + k is local wavenumber, k is modulation wavenumber, and 
a G C is called the growth rate. If Re (a) > 0, then the perturbed solution of the NLS 
equation grows exponentially. This is the criterion for the so-called Benjamin— Feir 
instability of a one-wave mode with modulation wavenumber n [3] . 



Substituting the function B(£,t) in (12) into (11) yields a pair of coupled equa- 

(13) 



tions that can be written in matrix form as follows 

ia-l3K 2 + 1 rl 7 rg W B x \ _ ( 



7^o -«7-/3K 2 + 7rg J \ B* 2 J \0 



Nontrivial solution to (13) can exist only if the determinant of the left hand side 



matrix is zero. This condition reads as follows 

a 2 = /?K 2 (2 7 rg -(3 k 2 ). (14) 

We have the following cases : 

• The growth rate a is real and positive if n 2 < 2 % rg . This corresponds to 
Benjamin— Feir instability. For specified values of k, the perturbation am- 
plitude is exponentially amplified in time [3] . 
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• The growth rate u is purely imaginary if k 2 > 2^Tq. This corresponds to a 
plane-wave solution that has bounded amplitude for all time [3]. 

Thus, the range of instability is given by 



< \k\ < \K crit 




It is easy to find that the 'strongest' instability occurs at K max = y jg r 0) where the 

maximum growth rate is cr max = 7 r o • Figure [l] shows the plot of growth rate a as 
a function of modulation wavenumber k for rg = /3 = 7 = 1. 



1.2 Growth rate plot for r = |3 = y = 1- 




0.5 1 1.5 

Modulation wavenumber K 

Figure 1: Growth rate a as a function of modulation wavenumber k for r$ = (3 = 
7 = 1. 

We can write the solution of the linearized NLS equation as 

^(e,r) = r e i ~ ir ° T [l + e^ T (B l e i ^ + B 2 e- lK % (16) 

with <j(k) = k a/2 (3 7 Tq — ft 2 k 2 . Since this solution is obtained from the linearized 
NLS equation, it is only valid if amplitudes are small. When the time is increased, 
the amplitude increases exponentially and the linearized theory becomes invalid. 
Therefore, we cannot use the solution from the linearized equation to investigate 
the behavior of the maximum amplitude increase due to the Benjamin-Feir instabil- 
ity. Fortunately, there exists an exact solution to the NLS equation that describes 
the exact behavior of the wave profile and that corresponds to the Benjamin-Feir 
instability. 

4 Modulational Instability 

In this section we investigate the relation between the maximum amplitude of a 
certain solution of the NLS equation and the modulation wavenumber k in the 
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instability interval. For simplicity, we choose the amplitude ro = 1, and the coeffi- 
cients f3 = 7 = 1. An exact solution, the so-called soliton on finite background, in 
short SFB, (sometimes also called the second most important solution), of the NLS 
equation is given by pQ : 



(k 2 — 1) cosh(<r(/c) r) + 



2-K 2 



cos(k£) + ia(n) sinh(cr(K) r) 



cos1i(<t(k) r) 



2-K 2 



cos(k^) 



(17) 



where < k < \[2 and ct(k) = — k 2 . 
For k = 1, we have 



, cos£ + iv^sinhr ir 
1>(S,t) := — e lT 



(18) 



y/2 cosh r — cos £ 

Figure [2] shows the plot of from (18) as a function of £ and r. Note that ip(£, r) 
is a 2-7T -periodic function with respect to £ variable. 



Solution M^,t)| for modulation waven urn bene - 1 




Figure 2: Plot of as a function of £ and r for modulation wavenumber k = 1. 



In the following, we analyze the solution (17) in detail. The behavior for this 
SFB as r — >■ ±00 is given by 



lim |V>(£, r)| = v 7 ^ 2 - l) 2 + (2 - k 2 ) = 1. 



(19) 



Because of this property the solution (17) is called as SFB. For a 'normal' soliton, 



the elevation vanishes at infinity : the 'normal' soliton is exponentially confined. For 
SFB, the solution is a similar elevation on top of the finite (nonzero) background 
level, here the normalized value 1. 
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Write the solution in the form ip^,r) = u(£,t) e lT , where u(£,t) where u de- 
scribes the amplitude and the exponential part expresses oscillations in time. Let us 
investigate the behavior of u(£, r) in time. For that consider the limiting behavior 
of d T u as r — > ±00. We have 



q u ia 2 (K) — <t{k) a/ 2 2 * 2 cos(k£) [k 2 sinh (<t(k) t) + ia(K) cosh (ct(k) t) 
d~r^ T) = = ^ 



If e^^, then 



and 



If £ = then 



and 



<9it 
97 



du 

du 
du 



cosh (<r(/c) r) — y^ 2 cos 

-2( K 2 +iCT(K))e~ ,J ( K ) T ifr»0, 
2(k 2 -k( K ))e ff(K)T ifr<0. 

a 4ia 2 (K)e- 2ff(K)T forr>0, 
« 4ia 2 (K)e 2<T(K)T forr<0. 



(20) 
(21) 

(22) 

(23) 
(24) 
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Maximum amplitude |\|/| at ^ = 0, X=0 
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Modulation wavenumber K 
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Figure 3: Plot amplitude maximum tp{^ = 0, r = 0; k) as a function of the modula- 
tion wavenumber k. 



Next, let us find the relation between the maximum amplitude of the exact 
solution (17) with modulation wavenumber k, where < k < y/2. The maximum 



8 



Natanael Karjanto, E. van Groesen, and Pearu Peterson 



value of the complex amplitude is at £ = (mod 2ir) and when r = 0. So, we have 



|^(0,0; K )\ 



1 



1 



i„2 



(25) 



Using the approximation \/l + a ~ 1 + for small a, and apply it to our formula 



(a = -k 2 ), we obtain 



lim |V>(0,0;/e)| = 3. 



As a result, the maximum factor of the amplitude amplification is 

|V(0,0;k)| 3 



(26) 



lim - 

k-+0 lim T ^ ±00 \ip{€,T; k)\ 1 



3. 



(27) 



Figure [3] shows the plot of the maximum amplitude |^| m ax as a function of modula- 
tion wavenumber k. 

To summarize, Figure [4] compares the plot of the dispersion relation f2 and its 
quadratic approximation (leads to the NLS equation), the growth rate a (which is 
related to Benjamin-Feir instability), and the maximum amplitude of as functions 
of wavenumber k. 




Wavenumber k 



Figure 4: Dispersion relation Q and its quadratic approximation, growth rate a 
(Benjamin-Feir instability), and maximum amplitude of ifi as functions of wavenum- 
ber k. The case corresponds to ko = 2ir. 
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5 Conclusion 

In this paper we modelled surface waves in wave tanks of hydrodynamics laboratories 
using the NLS equation. We analyzed the linearized NLS equation and obtained that 
its solutions have tendency to grow exponentially. We considered also the exact 
solution known as SFB of the NLS equation, that is the continuation of this linear 
instability. Using this, we found the maximum amplitude in space and time when 
the modulation wavenumber is in the interval of the Benjamin-Feir instability. The 
main result of this paper is that the maximum factor of the amplitude amplification 
due to the Benjamin-Feir instability is three. As we can see from Figure |4j the 
growth rate from the Benjamin-Feir instability does not determine the maximum 
amplitude amplification of the SFB. The results of this paper will be used in the 
further study of relating the Benjamin-Feir instability and the Phase- Amplitude 
equations. 
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